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Abstract 


A method for analysis of progressive failure in the Computational Structural 
Mechanics Testbed is presented in this report. The relationship employed in this anal- 
ysis describes the matrix crack damage and fiber fracture via kinematics-based 
volume-averaged variables. Damage accumulation during monotonic and cyclic 
loads is predicted by damage evolution laws for tensile load conditions . The imple- 
mentation of this damage model required the development of two testbed processors. 
While this report concentrates on the theory and usage of these processors , a com- 
plete list of all testbed processors and inputs that are required for this analysis are 
included. Sample calculations for laminates subjected to monotonic and cyclic loads 
were performed to illustrate the damage accumulation, stress redistribution , and 
changes to the global response that occur during the load history. Residual strength 
predictions made with this information compared favorably with experimental 
measurements. 


Introduction 

Laminated composite structures are susceptible to 
the development of microcracks during their operational 
lives. While these microcracks tend to aggregate in high 
stress regions and result in localized regions of reduced 
stiffness and strength, the microcracks can affect the glo- 
bal response of the structure. This change in the global 
structure in turn can create high stresses and increase 
damage accumulation in another part of the structure. 
Thus to accurately predict the structural response and 
residual strength of a laminated composite structure, the 
effects of the accumulating damage must be incorporated 
into the global analysis. The approach taken is to develop 
damage-dependent constitutive equations at the ply level. 
These equations are then employed in the development 
of the lamination equations from which the constitutive 
module of the structural analysis algorithm is con- 
structed. This algorithm is executed in a stepwise manner 
in which the damage-dependent ply-level results are used 
in the calculation of the global response for the next load 
step. This report will describe two Computational Struc- 
tural Mechanics (CSM) Testbed (COMET) processors 
that were developed for the performance of such an anal- 
ysis. A brief review of the theory behind the processors is 
first presented. The usage of these processors is then 
demonstrated. Since this analysis requires the use of 
other COMET processors, this report serves as a supple- 
ment to The Computational Structural Mechanics Test- 
bed User's Manual (ref. 1). 

It should be noted that the current damage model 
capability, computer code version 1.0, is limited to intra- 
ply matrix cracks and fiber fracture under tensile loads. 

Symbols and Abbreviations 

A laminate extensional stiffness matrix 

B laminate coupling stiffness matrix 


D 

laminate bending stiffness matrix 

dpara 

material parameter determined from 
experimental data 

E\\ 

lamina longitudinal modulus 

E 2 2 

lamina transverse modulus 

F 

applied force 

G 12 

lamina shear modulus 

k 

material parameter determined from 
experimental data 

N x 

applied load 

h 

material parameter determined from 
experimental data 

PDALC 

Progressive Damage Analysis of 
Laminated Composites 

R 

percent of maximum load 

'ply 

ply thickness 

u 

longitudinal extension 

u ° , v°, w° 

undamaged midplane displacements 

X , y 9 z 

displacement fields 

Vl2 

lamina Poisson’s ratio 


Damage-Dependent Constitutive Relationship 

The damage-dependent constitutive relationship 
employed in the COMET analysis is based on a contin- 
uum damage mechanics model proposed by Allen, 
Harris, and Groves (refs. 2 and 3). Rather than explicitly 
modeling each matrix crack in the material, the averaged 
kinematic effects of the matrix cracks in a representative 
volume are modeled by internal state variables. These 
internal state variables are defined by the volume- 
averaged dyadic product of the crack face displacement 



u\ and the crack face normal nj as proposed by 
Vakulenko and Kachanov (ref. 4): 

% = rjs'i"/* (l) 

where o!f is the second order tensor of internal state 
variables, V L is the local representative volume in the 
deformed state, and S is the crack surface area. This 
product can be interpreted as additional strains incurred 
by the material as a result of the internal damage. From 
micromechanics it has been found that the effects of 
matrix cracks can be introduced into the ply-level consti- 
tutive equation as follows (ref. 5): 

o L = [Q}{z l -o. m l } ( 2 ) 

where a L are the locally averaged components of stress, 
[Q] is the ply-level reduced stiffness matrix, and {e L } are 
the locally averaged components of strain. The laminate 
constitutive relationships are obtained by integrating the 
ply constitutive equations through the thickness of the 
laminate to produce 

W = [A]{e° L } + [B]{K L } + {f M } (3) 

{M} = [B]{e° l } + [D]{K t } + {g M } (4) 

where {TV} and {M} are the resultant force and moment 

vectors, respectively; [A], [2?], and [D] are the laminate 
extensional, coupling, and bending stiffness matrices, 
respectively (ref. 6); (e£) is the midplane strain vector; 
{k l } is the midplane curvature vector; and [f M } and 
{g M } are the damage resultant force and moment vectors 
for matrix cracking, respectively (ref. 7). The application 
of {f M } and \g M ) to the undamaged material will pro- 
duce midplane strain and curvature contributions equiva- 
lent to those resulting from the damage-induced 
compliance increase. 

As the matrix cracks accumulate in the composite, 
the corresponding internal state variables must evolve to 
reflect the new damage state. The rate of change of these 
internal state variables is governed by the damage evolu- 
tionary relationships. The damage state at any point in 
the load history is thus determined by integrating the 
damage evolutionary laws. Based on the observation that 
the accumulation of matrix cracks during cyclic loading 
is related to the strain energy release rate G in a power 
law manner (ref. 8), Lo, Allen, and Harris (ref. 9) have 
proposed the following evolutionary relationship for the 


internal state variable corresponding to the mode I (open- 
ing mode) of the matrix cracks: 

, M 

M - h 

da h 2 = ~if kGdN < 5 > 

The term da^ 22 reflects the changes in the internal 
state variable with respect to changes in the crack sur- 
faces. This term can be calculated analytically from a 
relationship that describes the average crack surface dis- 
placements in the pure opening mode (mode I) for a 
medium containing alternating 0° and 90° plies (ref. 5). 
The term G is the strain energy release rate calculated 
from the ply-level damage-dependent stresses. The mate- 
rial parameters k and h are phenomenological in nature 
and must be determined from experimental data (refs. 10 
and 11). Because k and h are assumed to be material 
parameters, the values determined from one laminate 
stacking sequence should be valid for other laminates as 
well. Since the interactions with the adjacent plies and 
damage sites are implicitly reflected in the calculation of 
the ply-level response through the laminate averaging 
process, equation (5) is not restricted to any particular 
laminate stacking sequence. 

When the material is subjected to quasi-static 
(monotonic) loads, the rate of change of the internal state 

variable is described by 

^22 

da M = { P^ £ 22 ~ e 22crit) if e 22 > e 22crit | 

^ e 22 <e 22criJ 

where £22crit * s critical tensile failure strain and (3 is a 
factor that describes the load carrying capability of the 
material after the critical tensile strain has been reached. 
Elastic perfectly plastic behavior is obtained by setting 
(3 = 1. A similar relationship is used to describe the ten- 
sile failure of the reinforcing fibers. The internal state 


variable for this mode of damage 
change is 

is o.¥ and its rate of 

l 22 

da L u = < 

\ Y^(£ n - £ ] lcrit ) 
1 

if E 11 >e llcrit| ^ 

1 1 

l o 

ife ll <e llcritJ 


where Icrit is the tensile fiber fracture strain and y is a 
factor describing the residual load carrying capability of 
the material after fiber fracture has occurred. 

Structural Analysis Formulation 

In order to simplify the formulation, it is expedient 
to consider the special case of symmetric laminates. With 
this case, the coupling stiffness matrix [B] becomes the 
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null matrix and the in-plane and out-of-plane laminate 
equations are decoupled. The laminate equations (3) and 
(4) are then substituted into the plate equilibrium equa- 
tions to yield the following governing differential equa- 
tions for the plate deformations: 
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These governing differential equations are integrated 
against variations in the displacement components 
to produce a weak form of the damage-dependent 
laminated-plate equilibrium equations. By substituting 
the corresponding displacement interpolation functions 
into the weak form of the plate equilibrium equations, the 
following equilibrium equations in matrix form are 
produced (ref. 12): 


Structural Analysis Scheme 

The nonlinear nature of the constitutive relationship 
and the progressive nature of the failure process requires 
that the analysis be performed in a stepwise manner as 
shown in figure 1 (from ref. 13). At each load step, the 
damage resultant forces and moments are determined for 
the current matrix and fiber damage state. The damage- 
induced force vector is then combined with the applied 
force vector. Nodal displacements are calculated with 
this combined force vector. The elemental stress result- 
ants are then determined. Finally, the ply-level stresses 
and strains are calculated as well as the damage evolution 
in each ply. This information is then used in the calcula- 
tions for the next load step. Because an iterative scheme 
to ensure equilibrium is not in place, each load step 
increment should be small enough to ensure an accurate 
solution. Since the effects of the matrix and fiber damage 
are represented as damage-induced force vectors, this 
formulation obviates the need to recalculate the elemen- 
tal stiffness matrices each time the damage state evolves. 
The fiber damage state is also used to determine the 
structural integrity of the component. Residual strength 
predictions can be made with this model by increasing 
the load or displacement at the boundary until fiber frac- 
ture is determined over a critical region of the compo- 
nent. This capability will be demonstrated in the 
following section entitled “Example Calculations.” 

The implementation of this analysis into the 
COMET code can be accomplished with the develop- 
ment of processors DRF (Damage Resultant Forces) and 
DGI (Damage Growth Increment). These processors, as 
with other COMET processors, are semi-independent 
computational modules that perform a specific set of 
tasks. Processor DRF first calculates the damage result- 
ant forces and moments and then incorporates them into 
the global force vectors. The second processor DGI post- 
processes the elemental stress resultants into ply-level 
stresses and strains by using the damage-dependent con- 
stitutive relationship. With this information, the proces- 
sor computes the damage evolution and updates the 
damage state for the next series of calculations. The 
remaining calculations can be performed with existing 
COMET processors. The following is a list in order of 
COMET processor executions for this analysis: 


[*]{S} = {F A } + {F M } (11) 

where [ K] is the element stiffness matrix, {5} is the dis- 
placement vector, {F^( is the applied force vector, and 
\Fm\ is the damage-induced force vector resulting from 
matrix cracking. Note that the effects of the internal dam- 
age now appear on the right side of the equilibrium equa- 
tions as damage-induced force vectors. 


1 . Procedure ES defines element parameters. 

2. Processor TAB defines joint locations, con- 
straints, and reference frames. 

3. Processor AUS builds tables of material proper- 
ties, section properties, and applied forces. 

4. Processor LAU forms constitutive matrix. 
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5. Processor ELD defines elements. 

6. Processor E initializes element data sets, and cre- 
ates element data sets. 

7. Procedure ES initializes element matrices. 

8. Procedure ES calculates element intrinsic stiff- 
ness matrices. 

9. Processor RSEQ resequences nodes for minimum 
total execution time. 

10. Processor TOPO forms maps to guide assembly 
and factorization of system matrices. 

1 1 . Processor K assembles system stiffness matrix. 

12. Processor INV factors system stiffness matrix. 

13. Continue. 

14. Processor DRF forms damage resultant force 
vectors. 

15. Processor SSOL solves for static displacements. 

16. Procedure STRESS calculates element stress 
resultants. 

17. Processor DGI calculates ply-level stresses and 
damage evolution. 

18. For next load cycle, go to step 13; else stop. 

The usage and theory behind each of the existing 
processors can be found in The Computational Structural 
Mechanics Testbed User's Manual (ref. 1). Information 
for processors DRF and DGI can be found in appendices 
A and B of this report, respectively. With the exception 
of processor DRF and DGI, other processors from the 
COMET processor library can be substituted into the list 
above to perform the tasks specified. 

Example Calculations 

Example calculations were conducted with COMET 
to illustrate the features of the progressive damage code. 
The first example demonstrates the effects of the evolv- 
ing matrix damage on a cross-ply laminated composite 
plate that was subjected to constant amplitude fatigue 
loads. The dimensions and boundary conditions for the 
laminated plate are shown in figure 2. This plate was dis- 
cretized into 24 four-node quadrilateral EX47 shell ele- 
ments (ref. 14). In this example, the plate has a [0/90] s 
laminate stacking sequence and the ply-level mechanical 
properties listed in table 1. These properties corre- 
sponded to those measured for IM7/5260 (ref. 11). A 
maximum load of 2500 lb/in at an tf-ratio of 0.1 was 
applied to the laminate. The COMET program and input, 
as well as a segment of the output, for this example can 


be found in the section entided “Progressive Failure 
Analysis Input” of appendix B. 

The predicted distribution of the mode I matrix crack 
damage in the 90° plies is shown in figure 3. The 
damage was greatest at the narrow end of the plate 
because the component of stress normal to the fiber was 
highest in this region. The higher stresses further trans- 
lated to a greater amount of energy available for the initi- 
ation and propagation of additional damage. This 
availability of energy was reflected in the damage evolu- 
tion along the length of the plate. However as damage 
accumulated in the plate, the stress gradient in the 90° 
plies became less steep (fig. 4). The similarity in stress 
resulted in relatively uniform changes to the damage 
state at higher load cycles. For this laminate stacking 
sequence, the load shed by the damaged 90° plies was 
absorbed by the 0° plies. The consequence of this load 
redistribution is an increase in the global displacements 
(fig. 5). The redistribution of load to the adjacent plies 
will affect the interlaminar shear stresses as well. This 
redistribution could create favorable conditions for the 
propagation of delamination. 

The second example examines the effects of damage 
accumulation during cyclic fatigue loads on the residual 
strength of notched laminates. For comparative purposes, 
unnotched laminates of similar dimensions were also 
examined. In this example, the notched laminates are ten- 
sion fatigue loaded for 100 000 cycles and then monoton- 
ically loaded to failure. The notched (central circular 
hole) laminate is shown in figure 6. Symmetry was 
assumed about the length and width of the laminate so 
that only a quarter of the laminate was modeled by the 
finite element model. This model, shown in figure 7, 
consisted of 153 four-node quadrilateral EX47 shell 
elements. Two laminate stacking sequences, a cross-ply 
[0/903] s a quasi -isotropic [0/±45/90]$, were consid- 
ered. These laminates possessed the same ply-level mate- 
rial properties as the first example. (See table 1.) The 
maximum fatigue loads employed in sample calculations 
are listed in table 2. The COMET program for the fatigue 
load portion of the calculation is similar to the example 
shown in appendix B. The residual strength portion of 
the program differed in that the monotonic matrix dam- 
age growth law, equation (6), is used in place of the 
fatigue damage growth law, equation (5). In addition, the 
applied load is incrementally increased with each load 
step to simulate a ramp up load input. Failure of the com- 
ponent is assumed to have occurred when the elements 
that have sustained fiber failure in the principal load car- 
rying plies span across the width of the laminate. The 
load at which this condition is satisfied is used to calcu- 
late the residual strength. At the present time, this struc- 
tural failure determination process is performed by the 
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analyst by tabulating the locations where fiber fracture is 
predicted. 

The COMET program where the laminate is first 
fatigue loaded then loaded monotonically to failure is 
listed in appendix C. In figure 8, the predicted stiffness 
loss for the open-hole geometry is compared to experi- 
mentally measured values of stiffness loss measured over 
a 4-in. gage length. The predicted residual strengths for 
the unnotched and open-hole geometries are shown along 
with experimental measurements in figure 9. 

The elastic perfectly plastic nonlinear behavior 
(y- 1) is a user specified assumption in the computer 
analysis. Other types of nonlinear materials behavior 
may also be selected by the user. For example, complete 
unloading (classical ply discount method y=0) can be 
assumed or any available strain softening law can be 
specified by the user. A comparison of the effect of the 
failure criteria on the longitudinal stresses in the 0° ply of 
the [0/±45/90] s laminate is shown in figure 10. Results 
for the undamaged stress state are compared to the redis- 
tribution loads (stresses) produced by the elastic per- 
fectly plastic criterion and the ply discount criterion at 
laminate failure. 

A systematic mesh refinement study was conducted 
for the quasi-isotropic laminate to determine if a numeri- 
cally converged analytical solution could be obtained. 
The analytical solutions for y = 1 converged after four 
successive refinements to the finite element mesh. The 
four meshes are shown in figure 1 1 and the numerical 
results of the convergence study are plotted in figure 12. 

Although this analysis considered only matrix crack- 
ing and fiber fracture, the results illustrate the effects of 
subcritical damage accumulation on the local and global 
response of a laminated composite. The inclusion of 
other damage modes such as delamination and compres- 
sion failure mechanisms will provide a more complete 
picture of the failure process. Since matrix cracking usu- 
ally precedes these two modes of damage, the present 
analysis can be employed to determine the initiation and 
propagation of these other modes of damage. Finally, the 
introduction of failure criteria for additional modes of 
damage would enable the prediction of the progressive 


failure process up to catastrophic failure of laminated 
composite structures (ref. 14). 

Concluding Remarks 

This report describes a progressive failure analysis 
for laminated composites that can be performed with the 
Computational Structural Mechanics (CSM) Testbed 
(COMET) finite element code. The present analysis uses 
a constitutive model that describes the kinematics of the 
matrix cracks via volume averaged internal state vari- 
ables. The evolution of these internal state variables is 
governed by an experimentally based damage evolution- 
ary relationship. The nonlinearity of the constitutive 
relationship and of the damage accumulation process 
requires that this analysis be performed incrementally 
and possibly iteratively. 

Two processors were developed to perform the nec- 
essary calculations associated with this constitutive 
model. In the analysis scheme, these processors were 
called upon to interact with existing COMET processors 
to perform the progressive failure analysis. This report, 
which serves as a guide for performing progressive fail- 
ure analysis on COMET, provides a brief background on 
the constitutive model and the analysis methodology in 
COMET. The description and usage of the two progres- 
sive failure processors can be found in the appendices of 
this report. These appendices are meant to supplement 
the COMET User’s Manual. 

The results from the example problems illustrate the 
stress redistribution that occurred during the accumula- 
tion of matrix cracks and fiber fracture. This stress redis- 
tribution in turn influenced the damage evolution 
characteristics, the global displacements, and the residual 
strengths. It should be noted that the current damage 
model capability, computer code version 1.0, is limited 
to intraply matrix cracking and fiber fracture under ten- 
sile load conditions. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
July 24, 1996 


5 



Appendix A 


Processor DRF 

Al. General Description 

This processor calculates the damage resultant forces and moments caused by matrix cracking in laminated compos- 
ites. These resultant forces and moments when applied to an undamaged laminate will produce an equivalent amount of 
displacements and curvatures to those resulting from the matrix crack surface kinematics in a damaged laminate. This 
enables an analysis of the response of a damaged laminate without having to update the stiffness matrix each time the 
damage state changes. Matrix crack damage is modeled in this processor by volume averaged crack surface kinematics 
that use internal state variables (refs. 2 and 3). 

Processor DRF and processor DGI, which is described in appendix B, were developed to perform progressive fail- 
ure analysis of quasi-static and fatigue loaded laminates in the Computational Structural Mechanics (CSM) Testbed 
(COMET). Analyses from these processors are stored in two formats. One is in standard format that is accessed by open- 
ing the output file. The other is a data set, which is stored in a testbed data library, and provides data to processors and 
post-processors (ref. 1). In this analysis, processor DRF is used in conjunction with COMET analysis processors to 
determine the static displacement and elemental stress resultants for a laminated composite structure containing matrix 
crack damage. Processor DGI then calculates the damage-dependent ply stresses. The damage state is updated based on 
the ply stresses and this procedure is repeated for the next load cycle. 

Al.l. Damage-Dependent Constitutive Relationship 

In this processor, the effects of the matrix cracks are introduced into the ply-level constitutive equations as follows 
(ref. 3): 


{° L } = [(2]{e L -«f} (Al) 

where {a^} are the locally averaged components of stress, [Q] is the ply-level reduced stiffness matrix, and {e L } are 
the locally averaged components of strain. The variables {a^} are the components of the strain-like internal state 
variable for matrix cracking and are defined by 


a 


M 

h 


J 

V 


L 


\ s u ‘ n j ds 


(A2) 


where V L is the volume of an arbitrarily chosen representative volume of ply thickness that is sufficiently large that 
do not depend on V L , is the crack opening displacement, rij is the component of the vector normal to the crack face, 
and S is the surface value of the volume V L . The present form of the model assumes that a^, the internal state variable 
representing the mode I matrix crack opening, is the only nonzero component. 


A1.2. Damage-Dependent Laminate Equations 

The ply-level strains are defined as follows: 


^ Lxx ^Lxx zK Lxx 

(A3) 

Z Lyy = £ Lyy ~ zK Lyy 

(A4) 

£ Lxy = £ Lxy~ K Lxy 

(A5) 
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where e° L and k l are the midplane strains and curvatures, respectively. The aforementioned ply strains are then substi- 
tuted into equation (Al) to produce the ply-level stresses. Damage-dependent lamination equations are obtained by 
integrating these ply stresses through the thickness of the laminate (ref. 15). Next, the stiffness matrix in the laminate 
equation is inverted to produce 


A B 
B D 


~ U N-f 


M-g 


M 


(A6) 


where [A], [5], and [D] are, respectively, the undamaged laminate extensional, coupling, and bending stiffness matrices. 
They are defined by the following equations from reference 6: 


[A] = Y^[QUz k -z kA ) (A7) 

*=1 

n 

[*] = \^lQh(z 2 k -z 2 k . ,) (A8) 

k= 1 
n 

[D] = l£[CWz5-zL) (A9) 

k=\ 

where [Q] k is the transformed reduced elastic modulus matrix for the kt h ply in laminate coordinates. In equation (A6), 
N is the component of the resultant force per unit length and M is the component of the resultant moment per unit length. 
The variables { f M } and { g M } represent the contribution to the resultant forces and moments from matrix cracking and 
are calculated from 


{/*} = -X[G]*(z t -z*_ 1 ){«£}* (A10) 

*=1 

is M } = -l^mizl-zl i ){ <xf } * (aid 

k = 1 

where {a**}* contains the matrix cracking internal state variables for the k\h ply. Thus given the forces N and moments 
A/, as well as the damage variables in each ply, equation ( A6) can be used to calculate the midsurface strains z L and cur- 
vature 

A2. Processor Syntax 

This processor uses keywords and qualifiers along with the CLIP command syntax (ref. 1). Two keywords are rec- 
ognized: SELECT and STOP. 
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A2.1. Keyword SELECT 


This keyword uses the qualifiers listed below to control the processor execution. 


Qualifier 

Default 

Meaning 

LIBRARY 

1 

Input and output library. 

ELEMENT 

ALL 

Element type (EX47, EX97) used in the analysis. Default is all 
element types found in LIBRARY. 

SREF 

1 

Stress reference frame. Stress resultants may have been computed in 
the element stress/strain reference frame (SREF = 0) or in one of 
three alternate reference frames. For SREF = 1, the stress/strain 
x-direction is coincident with the global ^-direction. For SREF = 2, 
the stress/strain ^-direction is coincident with the global y-direction. 
For SREF = 3 the stress/strain ^-direction is coincident with the 
global z-direction. Note that the processor currently must have the 
stress/strain coincident with the global ^-direction fSREF = 1 h 

PRINT 

1 

Print flag. May be 0, 1, or 2; 2 results in the most output. 

MEMORY 

2000 000 

Maximum number of words to be allocated in blank common. This 
is an artificial cap on memory put in place so that the dynamic 
memory manager does not attempt to use all of the space available 
on the machine in use. 

DSTATUS 

1 

Damage state flag. If no damage, DSTATUS = 0. If matrix cracking 
(cyclic load), DSTATUS = 1. If matrix cracking (monotonic load), 
DSTATUS = 22222. 

X FACTOR 

0.0 

Increases the specified applied forces by this factor at every load 
step. This qualifier is used in the residual strength calculations. 


A2.2. Keyword STOP 

This keyword has no qualifiers. 

A3. Subprocessors and Commands 

Processor DRF does not have subprocessors. 
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A4. Processor Data Interface 


A4.1. Processor Input Data Sets 

Several data sets, which are listed below, are used as input for processor DRF. 


Input data set 

Contents 

ELTS.NAME 

Element names 

OMB.DATA.1.1 

Material properties including strain allowables 

LAM.OMB.*.* 

Laminate stacking sequence 

ES. SUMMARY 

Various element information 

PROP.BTAB.2.102 

ABD matrix 

WALL.PROP.1.1 

Shell wall data set 

DIR.xxxx.*.* 

Element directory data set 

DEF.xxxx.*.* 

Element definition (connectivity) data set 

ISV.xxxx.*.* 

Internal state variable data set 

xxxx.EFIL.*.* 

Element nodal coordinates and transformations 

APPL.FORC 

Applied force and moments at joints 


A4.2. Processor Output Data Sets 

These data sets are used as output for processor DRF. 


Output data set 

Contents 

APPL.FORC 

Applied force and moments at joints 

DFCT.xxxx.*.* 

Temporary damage resultant force data set 

DRFC.xxxx.*.* 

Damage resultant force data set 


A5. Limitations 

Only EX47 and EX97 elements implemented with the generic element processor ESI will be processed by proces- 
sor DRF. All other elements will be ignored. The stress reference frame must be coincident with the global ^-direction. 

A6. Error Messages 

Fatal errors will occur when any of the required data sets are missing from the input data library or when the stress 
resultants at the integration points are missing. (See section A4.1.) 

Warning messages will be written and execution will continue when there is a missing or unreadable keyword or 
qualifier or if any of the original SPAR elements are encountered. 

A7. Usage Guidelines and Examples 

A7.1. Program Organization 

The following list illustrates the organization of a progressive failure analysis that uses COMET. Because of the 
nonlinear nature of the damage-dependent constitutive equation, this analysis is performed in a stepwise manner. With 
the exception of processors DRF and DGI, all COMET processors can be employed to perform the specified tasks. 
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1 . Procedure ES defines element parameters 

2. Processor TAB defines joint locations, constraints, and reference frames 

3. Processor AUS builds tables of material and section properties and applied forces 

4. Processor LAU forms constitutive matrix 

5. Processor ELD defines elements 

6. Processor E initializes and creates element data sets 

7. Procedure ES initializes element matrices 

8. Procedure ES calculates element intrinsic stiffness matrices 

9. Processor RSEQ resequences nodes for minimum total execution time 

10. Processor TOPO forms maps to guide assembly and factorization of system matrices 

1 1 . Processor K assembles system stiffness matrix 

12. Processor INV factors system stiffness matrix 

13. Continue 

14. Processor DRF forms damage resultant force vectors 

15. Processor SSOL solves for static displacements 

16. Procedure STRESS calculates element stress resultants 

17. Processor DGI calculates ply-level stresses and damage evolution 

18. For next load cycle, go to step 13; else stop 

A7.2. Progressive Failure Analysis Input and Output 

Please refer to processor DGI in appendix B for an example. 

A8. Structure of Data Sets Unique to Processor DRF 

A8.1. DRFC.xxxx 

This data set is created by processor DRF and uses the SYSVEC format. See APPL.FORC.iset. 1 . This data set con- 
tains the damage resultant forces and moments corresponding to the given matrix cracking damage state. 

A8.2. DFCT.xxxx 

Data set DFCT.xxxx is created by processor DRF and uses the SYSVEC format. See APPL.FORC.iset. 1 . This data 
set contains the damage resultant forces and moments from the previous load step and is used to restore the applied force 
vector to the initial value. 

A8.3. ISV.xxxx 

This data set contains the matrix cracking internal state variables at each layer. The xxxx is the element name. The 
data are stored in a record named ALP AM. 1 . This record contains n items, where 

n = nlayer * nintgpt * nelt 

and nlayer is the number of layers in the model, nintgpt is the number of integration points for element, and nelt 
is the number of elements. 
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The data are stored in the following order: 


M 

1. a, is internal state variable associated with fiber fracture. 
L \\ 


2. a/ is internal state variable associated with mode I opening of the matrix crack. 
l 22 


3. a, is internal state variable associated with mode II opening of the matrix crack. 
L \2 


The data storage occurs for every layer, every integration point, and every element. 
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Appendix B 


Processor DGI 

Bl. General Description 

Processor DGI predicts the evolution of matrix crack damage in laminated composites for monotonical loads and 
cyclic fatigue loads. The processor also calculates fiber fracture under tensile load conditions. The matrix crack damage 
is represented in this processor by volume-averaged crack surface kinematics that use internal state variables (refs. 2 
and 3). The evolution of these internal state variables is governed by a phenomenological growth law. 

This processor was designed to perform progressive failure analysis of laminated composite structures in the Com- 
putational Structural Mechanics (CSM) Testbed (COMET). At each load cycle, the elemental stress resultants for a lam- 
inated composite structure are obtained from COMET with the effects of matrix crack damage accounted for by 
processor DRF. Processor DGI then postprocesses this information and uses the ply-level stresses to determine the evo- 
lution of matrix crack damage in each ply of the laminate. This procedure is repeated until the specified number of load 
cycles has been reached. 

Bl.l. Damage-Dependent Constitutive Relationship 

In this processor, the effects of the matrix cracks are introduced into the ply-level constitutive equations as follows 
(ref. 5): 

{a L } = [Q]{e L - a“] (Bl) 

where {a L } is the locally averaged component of stress, [Q] is the ply-level reduced stiffness matrix, and {e L } are the 
locally averaged components of strain. The { a } are the components of the strain-like internal state variable for matrix 
cracking and are defined by 


% = \ rJ s u i n j dS (B2 > 

where V L is the volume of an arbitrarily chosen representative volume of ply thickness that is sufficiently large that 
do not depend on V L , w; are the crack opening displacements, and nj are the components of the vector normal to the 
crack face. The present form of the model assumes that a^ 2 , the internal state variable representing the mode I matrix 
crack opening, is the only nonzero component. 

For a uniaxially loaded medium containing alternating 0° and 90° plies, has been found from a micromechan- 
ics solution to be related to the far field normal force and crack spacing as follows (ref. 5): 



£ 

2r 


4 

64 ^ 2222 


(B3) 


where 


5 = II 2 1 73 4 (B4) 

m=ln=l C 22222 ( < 2m ~ *) ( 2/? “ 1 ) + C nn( a/t ) ( 2n ~ J ) 

p is the force per unit length that is applied normal to the fibers and It and 2 a are the layer thickness and crack spacing, 
respectively. The C 2222 the modulus in the direction transverse to the fibers and the C i212 the in-plane shear modu- 
lus. Both moduli are the undamaged properties. 


12 



B1.2. Damage Evolution Relationship 


Equation (B3) is used when the matrix crack spacing is known in each ply of the laminate. Since it is usually neces- 
sary to predict the damage accumulation and response for a given load history, damage evolutionary relationships must 
be utilized to determine the values of the internal state variables. The following relationship was used for the rate of 
change of the internal state variable a ^ in each ply during fatigue loading conditions (ref. 9): 

da M 

da? = —?^kG n dN (B5) 

dS 

where do!^ 22 describes the change in the internal state variable for a given change in the crack surface areas, k and h 
are material parameters (refs. 10 and 1 1), N is the number of load cycles, and G is the damage-dependent strain energy 
release rate for the ply of interest and is calculated from the following equation: 


G = 


V L C ijkl\ 





f a Lkl 

l dS 


(B6) 


where V L is the local volume. Interactions with the adjacent plies will result in ply strains , which are affected by the 
strains in adjacent plies. Thus, the strain energy release rate G in each ply will be implicitly reflected in the calculation of 
the ply-level response, so that equation (B5) is not restricted to a particular laminate stacking sequence. Substituting 
equation (B6) in equation (B5) and integrating the result in each ply over time gives the current damage state in each ply 
for any fatigue load history. 

When the material is subjected to monotonically increasing loads, the rate of change of the internal state variable 
a, is described by 



P^( £ 22 " £ 22crit) £ 22 > E 22crit j 

0 if £22 <e 22criJ 


(B7) 


where e 2 2crit * s critical tensile failure strain and p is a factor that describes the load carrying capability of the material 
after the critical tensile strain has been reached. Elastic perfectly plastic behavior is obtained by setting P = 1. A similar 
relationship is used to describe the tensile failure of the reinforcing fibers. The internal state variable for this mode of 
damage is and its rate of change is 



^( £ ll- £ llcrit 
0 


it) ife ii> e 


1 lent 


if e n <e Ucrit 


(B8) 


where £j j cr j t is the tensile fiber fracture strain and y is a factor describing the residual load carrying capability of the 
material after fiber fracture has occurred. 


B2. Processor Syntax 

This processor uses keywords and qualifiers along with the CLIP command syntax. Two keywords are recognized: 
SELECT and STOP. 
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B2.1. Keyword SELECT 

This keyword uses the qualifiers listed below to control the processor execution. 


Qualifier 

Default 

Meaning 

LIBRARY 

i 

Input and output library. 

ELEMENT 

ALL 

Element type (EX47, EX97) used in the analysis. Default is all element 
types found in LIBRARY. 

LOADJSET 

1 

Load set; i of input data set STRS .xxxx.i.j. 

SREF 

1 

Stress reference frame. Stress resultants may have been computed in 
the element stress/strain reference frame (SREF = 0) or in one of three 
alternate reference frames. For SREF = 1, the stress/strain x-direction is 
coincident with the global jc-direction. For SREF = 2 , the stress/strain 
^-direction is coincident with the global y-direction. For SREF = 3 the 
stress/strain ^-direction is coincident with the global z-direction. Note that 
the processor currently must have the stress/strain coincident with the 
global ^-direction fSREF = 11. 

PRINT 

1 

Print flag. May be 0, 1, or 2; 2 results in the most output. 

STEP 

0 

Step number in nonlinear analysis (i.e., i in the STRS.xxxx.i.O data set for 
nonlinear analysis). 

MEMORY 

2 000 000 

Maximum number of words to be allocated in blank common. This is an 
artificial cap on memory put in place so that the dynamic memory manager 
does not attempt to use all of the space available on the machine in use. 

DSTATUS 

1 

Damage state flag. If no damage, DSTATUS = 0. If matrix cracking (cyclic 
load), DSTATUS = 1 . If matrix cracking (monotonic load), DSTATUS = 
22222. 

INCSIZE 

1.0 

Increment size used in damage growth law. 

NCYCLE 

1 

Cycle number. 

NINCR 

1 

Increment number. 


B2.2. Keyword STOP 

This keyword has no qualifiers. 

B3. Subprocessor and Commands 

None. Processor DGI does not have subprocessors. 
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B4. Processor Data Interface 


B4.1. Processor Input Data Sets 

Several data sets, which are listed below, are used as input for processor DGI. 


Input data set 

Contents 

ELTS.NAME 

Element names 

STR S.xxxxJJ 

Element stress resultants. Record named INTEG_PTS must exist. 

OMB. DATA. 1.1 

Material properties including strain allowables 

LAM.OMB.*.* 

Laminate stacking sequence 

ES. SUMMARY 

Various element information 

PROP.BTAB.2. 1 02 

ABD matrix 

ISV.xxxx.*.* 

Internal state variable data set 

DEF.xxxx.*.* 

Element definition (connectivity) data set 

WALL. PROP. 1 . 1 

Shell wall data set 

DIR.xxxx.*.* 

Element directory data set 

DGP.DATA. 1 . 1 

Damage growth law parameters data set 


B4.2. Processor Output Data Sets 


Output data set 

Contents 

ISV.xxxx.*.* 

Internal state variable data set 

PDAT.xxxx 

Ply-level stresses, strains, and damage state 


B5. Limitations 

Only EX47 and EX97 elements implemented with the generic element processor ES 1 will be processed by proces- 
sor DGI. All other elements will be ignored. Currently, the stress reference frame must be coincident with the global 
x-direction. 

B6. Error Messages 

Fatal errors will occur when any of the required data sets are missing from the input data library or when integration 
point values of the stress resultants are missing. (See section B4.1.) 

Warning messages will be written and execution will continue when there is a missing or unreadable keyword or 
qualifier or when any of the original SPAR elements are encountered. 
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B7. Usage Guidelines and Examples 


B7.1. Organization of Progressive Failure Analysis on Testbed 

The organization of the COMET processors for a progressive failure analysis is shown below. The nonlinear nature 
of the damage-dependent constitutive equation requires that this analysis be performed in a stepwise manner. With the 
exception of processors DRF and DGI, any COMET processors can be called upon to perform the required tasks. 

1 . Procedure ES defines element parameters 

2. Processor TAB defines joint locations, constraints, reference frames 

3. Processor AUS builds tables of material and section properties and applied forces 

4. Processor LAU forms constitutive matrix 

5. Processor ELD defines elements 

6. Processor E initializes and creates element data sets 

7. Procedure ES initializes element matrices 

8. Procedure ES calculates element intrinsic stiffness matrices 

9. Processor RSEQ resequences nodes for minimum total execution time 

10. Processor TOPO forms maps to guide assembly and factorization of system matrices 

1 1 . Processor K assembles system stiffness matrix 

12. Processor INV factors system stiffness matrix 

13. Continue 

14. Processor DRF forms damage resultant force vectors 

15. Processor SSOL solves for static displacements 

16. Procedure STRESS calculates element stress resultants 

17. Processor DGI calculates ply-level stresses and damage evolution 

18. For next load cycle, go to step 13; else stop 

B7.2. Progressive Failure Analysis Input 

The following list illustrates the input from a progressive failure analysis. The uniaxially tensile-loaded tapered lam- 
inated plate, which is described in the main body of this report, is being solved (fig. 2). The list contains the main pro- 
gram plus a procedure file to perform the calculations for each load cycle. 

# 

#@\$-me 

# 

cp $CSM__PRC/proclib . gal proclib.gal .Copy procedure library 
chmod u+w proclib.gal 
testbed << \endinput 

*set echo=off .Do not echo input 

*set plib=28 

*open 28 proclib.gal /old .Open procedure library 

*open/new 1 qoutput.101 .Open output library 

tapered panel 

EX47 4 node quad elements 

24 nodes, 14 elements 
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*add pffc.clp .Add procedure for repeating calculations 

*def/a es_name = EX47 .Element name 

*def/a es_proc = ESI .Element processor name 

*call ES ( function= ' DEFINE ELEMENTS’; es_proc = <es_proc> 
es_name=<es_name> ) 


[ xqt TAB 

START 24 .24 nodes 

JOINT LOCATIONS .Enter joint locations 

1 0.0 0.0 0.0 20.0 2.0 0.0 813 

8 0.0 10.0 0.0 20.0 8.0 0.0 

CONSTRAINT DEFINITION 1 .Constraints: 


zero 1 , 2 , 3 , 4 , 5 : 1,17,8 
zero 6: 1,24 


.Fixed end 

.Suppress drilling DOF 


[xqt AUS 

SYSVEC : appl fore 


.Create input datasets 
.Applied forces 


1=1 : J=8 : 3750.0 

1=1 : J=16 : 7500.0 

1=1 : J=24 : 3750.0 

TABLE (NI=16,NJ=1) : OMB DATA 1 1 .Ply-level material data 
IM7/5260 
1=1, 2,3, 4, 5 

J=1 : 22 . 162E+6 0.333 1.262E+6 0.754E+6 0.754E+6 
1=6,7, 8, 9 

J=l: 0.754E+6 1.0E-4 1.0E-4 0.01 
1=10, 11, 12,13,14,15, 16 

J=l: 0.0 0.0 0.0 0.014 0.014 0.014 0.0 

TABLE (NI=3 ,NJ=3 , ITYPE=0) : LAM OMB 1 1 .Section properties 


J = 1 : 

1 

0.006 

0 . 0 

J=2 : 

1 

0.012 

90.0 

J = 3 : 

1 

0.006 

0.0 


TABLE (NI=3 , NJ=1 , ITYPE=0 ) : DGP DATA 1 1 .Damage evolution data 
J=1 : 1.1695 5.5109 3.8686E-7 


[xqt LAU 
ONLINE=2 
[xqt ELD 

<es_expe_cmd> 

NSECT = 1 : SREF=1 
1 
2 

3 

4 

5 

6 
7 
9 


2 

3 

4 

5 

6 

7 

8 

10 


.Create constitutive matrix 
. Define connectivity 


10 

9 

11 

10 

12 

11 

13 

12 

14 

13 

15 

14 

16 

15 

18 

17 


17 



datasets 


10 

11 

11 

12 

12 

13 

13 

14 

14 

15 

15 

16 


[xqt E 
stop 

*call ES ( function^ ' INITIALIZE 1 ) 

*call ES (function^ ’ FORM STIFFNESS /MATL ’ ) 
[xqt RSEQ 

reset maxcon=12 
[xqt TOPO 
[xqt K 
[xqt INV 
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18 

20 

19 

21 

20 

22 

21 

23 

22 

24 

23 


. Initialize element 

.Initialize element matrices 
.Form stiffness matrices 
. Resequence 

. Create maps 

.Assemble global stiffness matrix 
.Factor the global stiffness matrix 


*def/i ns_overwrite=<true> 

Call procedure to perform calculations at each cycle 

*call PFFC ( es_proc=<es_jproc> ; es_name=<es_name> ; 

N_f cycl=l ; N_lcycl=2000 ; N_cylinc=5 ; -- 
NPRT=1 0 0 ) 


*pack 1 
[xqt exit 
\endinput 


B7.2.1. Procedure to perform loop through calculations for each load cycle (file name pffc.clp) 

^procedure PFFC ( es_proc ; es_name ; 

N_f cycl ; N_lcycl ; N_cylinc ; -- 
NPRT ) 


N_fcycl : first fatigue cycle 
N_lcycl : last fatigue cycle 
N_cylinc: cycle increment 
NPRT: output storage cycle increment 

begin loop here 

*set echo=on,ma 
*set echo=off 

*def icount = 0 .Initialize print counter 

*D0 : CYCLOOP $NCYL = <[N_fcycl]>, <[N_lcycl]> / <[N_cylinc]> 
*def icount = ( <icount> + 1 ) 

*if < <icount> /eq <[NPRT]> /or <$NCYL> /eq 1 > /then 
*def iprint = 1 
*def icount = 0 
*else 
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0 


*def iprint = 

*endif 

*def delinc = < [N__cylinc ] > 

[xqt DRF .Calculate damage resultant forces 

select /PRINT = 0 
stop 

[xqt SSOL .Solve for static displacements 


Calculate elemental stress resultants 
*call STRESS (directional; location^ INTEG_PTS; print=<f alse> ) 
*if < < I PRINT> /eq 1 > /then 

[xqt VPRT .Print static displacements 

format = 4 
print STAT DISP 
stop 

[xqt DGI 

select /PRINT = 2 
select /INC_SIZE = <delinc> 
select /N_CYCLE = <$NCYL> 
select /NINCR = 1 
select /NINCR = <$SNCYL> 
stop 

*endif 

*if < <IPRINT> /ne 1 > /then 
[xqt DGI 

select /PRINT = 0 
select /INC_SIZE = <delinc> 
select /N_CYCLE = <$NCYL> 
select /NINCR = 1 
select /NINCR = <$SNCYL> 
stop 


.Calculate ply-level stresses, 
.strains, and damage evolution 


.Calculate ply-level stresses, 
.strains, and damage evolution 


*endi f 

*if < < IPRINT> /eq 1 > /then .Store datasets for post processing 

*copy 1, PLYDT . <ES_NAME> . <$NCYL> . 1 = 1, PDAT . <ES_NAME> . * 

*copy 1, DISP . <ES_NAME> . <$NCYL> . 1 = 1, STAT. DISP.* 

*endif 


CYCLOOP 

*set echo=off 
end 
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B7J. Progressive Failure Analysis Output 

The following is a partial list of a progressive failure analysis output produced by processor DGI. Data for post- 
processing are stored in data set PLYDT.xxxx.xxx.l . 

** BEGIN DGI ** DATA SPACE= 2000000 WORDS 
CYCLE NUM. = 496 


ELEMENT NUMBER 1 TYPE EX47 

EVALUATION (INTG) POINT NUMBER 1 

REFERENCE SURFACE STRAINS AND CURVATURES 

EO-X E0-Y E0-XY K-X K-Y K-XY 

0 . 4619E-02 -0 . 6946E-04 0.1180E-02 O.OOOOE+OO O.OOOOE+OO 0.0000E+00 

COMBINED BENDING AND MEMBRANE STRESSES, STRAINS, AND MATRIX CRACK DAMAGE 
VARIABLE 

FOR EACH LAYER OF ELEMENT 1 TYPE EX47 


LAYER 

THETA 

SIG-l 

SIG-2 

TAU-12 

STRAIN-1 

STRAIN-2 

GAMMA-12 

1 

0 . 

0. 103E+06 

0 . 187E+04 

0 . 89 0E+03 

0.462E-02 

-0.695E-04 

0 . 118E-02 

2 

90 . 

0 . 384E+03 

0 „ 57 8E+04 

-0.890E+03 

-0 . 695E- 04 

0 . 462E-02 

-0. 118E-02 

3 

0 . 

0 . 103E+06 

0 . 187E+04 

0.89 0E+03 

0 . 462E-02 

-0.695E-04 

0 . 118E-02 


LAYER ALPM-11 ALPM-22 ALPM-12 

1 0 . 000E+00 0 . 122E-11 O.OOOE+OO 

2 O.OOOE+OO 0 . 473E-04 O.OOOE+OO 

3 O.OOOE+OO 0 . 122E-11 O.OOOE+OO 


** BEGIN DGI ** DATA SPACE= 2000000 WORDS 
CYCLE NUM. = 996 


ELEMENT NUMBER 1 TYPE EX47 

EVALUATION (INTG) POINT NUMBER 1 

REFERENCE SURFACE STRAINS AND CURVATURES 

EO-X E0-Y E0-XY K-X K-Y K-XY 

0 . 4623E-02 -0.6882E-04 0.1183E-02 0.0000E+00 0.0000E+00 O.OOOOE+OO 

COMBINED BENDING AND MEMBRANE STRESSES, STRAINS, AND MATRIX CRACK DAMAGE 
VARIABLE 

FOR EACH LAYER OF ELEMENT 1 TYPE EX47 


LAYER 

THETA 

SIG-l 

SIG-2 

TAU-12 

STRAIN-1 

STRAIN-2 

GAMMA-12 

1 

0 . 

0 . 103E+06 

0 . 187E+04 

0. 892E+03 

0 .462E-02 

-0. 688E-04 

0.118E-02 

2 

90. 

0 .382E+03 

0 . 573E+04 

-0. 892E+03 

-0.688E-04 

0.462E-02 

-0 . 118E-02 

3 

0 . 

0 . 103E+06 

0 . 187E+04 

0 . 892E+03 

0.462E-02 

- 0 . 688E-04 

0.118E-02 


LAYER ALPM-11 ALPM-22 ALPM-12 

1 O.OOOE+OO 0 . 246E-11 O.OOOE+OO 

2 O.OOOE+OO 0.901E-04 0.000E+00 

3 0.000E+00 0.246E-11 0.000E+00 


** BEGIN DGI ** DATA SPACE= 2000000 WORDS 
CYCLE NUM. = 1496 


ELEMENT NUMBER 1 TYPE EX47 

EVALUATION (INTG) POINT NUMBER 1 
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REFERENCE SURFACE STRAINS AND CURVATURES 

EO-X EO-Y EO-XY K-X K-Y K-XY 

0 . 4625E-02 - 0 . 6839E- 04 0.1184E-02 O.OOOOE+OO O.OOOOE+OO O.OOOOE+OO 

COMBINED BENDING AND MEMBRANE STRESSES, STRAINS, AND MATRIX CRACK DAMAGE 
VARIABLE 

FOR EACH LAYER OF ELEMENT 1 TYPE EX47 


LAYER 

THETA 

SIG-l 

SIG-2 

TAU-12 

STRAIN-1 

STRAIN-2 

GAMMA-12 

1 

0 . 

0 . 103E+06 

0 . 187E+04 

0 . 89 3E + 03 

0 .463E-02 

-0 . 684E-04 

0 . 118E-02 

2 

90 . 

0 . 376E+03 

0 . 568E+04 

-0 . 893E+03 

-0 . 684E-04 

0 . 463E-02 

-0 .118E-02 

3 

0 . 

0 . 103E+06 

0 . 187E+04 

0.893E+03 

0 .463E-02 

-0 . 684E-04 

0.118E-02 


LAYER ALPM-11 ALPM-22 ALPM-12 

1 0 . OOOE+OO 0 . 372E-11 0 . OOOE+OO 

2 0. OOOE+OO 0 . 129E-03 0. OOOE+OO 

3 0. OOOE+OO 0 . 372E-11 0. 000E+00 


** BEGIN DGI ** DATA SPACE= 2000000 WORDS 
CYCLE NUM. = 1996 


ELEMENT NUMBER 1 TYPE EX47 

EVALUATION ( INTG) POINT NUMBER 1 

REFERENCE SURFACE STRAINS AND CURVATURES 

EO-X E0-Y EO-XY K-X K-Y K-XY 

0 . 4627E-02 -0.6806E-04 0.1185E-02 O.OOOOE+OO 0.0000E+00 O.OOOOE+OO 

COMBINED BENDING AND MEMBRANE STRESSES, STRAINS, AND MATRIX CRACK DAMAGE 
VARIABLE 


FOR 

EACH 

LAYER OF ELEMENT 1 

TYPE EX47 




LAYER 

THETA 

SIG-l 

SIG-2 

TAU-12 

STRAIN-1 

STRAIN-2 

GAMMA-12 

1 

0 . 

0 . 103E+06 

0 . 187E+04 

0 . 894E+03 

0.463E-02 

-0.681E-04 

0 . 119E-02 

2 

90. 

0.37 0E+03 

0 . 564E+04 

-0 .894E+03 

-0 .681E-04 

0 .463E-02 

-0 .119E-02 

3 

0 . 

0 . 103E+06 

0 . 187E+04 

0.894E+03 

0.463E-02 

-0 . 681E-04 

0 . 119E-02 


LAYER ALPM-11 ALPM-22 ALPM-12 

1 0. 000E+00 0 . 500E-11 0. 000E+00 

2 0. 000E+00 0 . 164E-03 0. OOOE+OO 

3 0. 000E+00 0 . 500E-11 0. OOOE+OO 


B8. Structure of Data Sets Unique to Processor DGI 
B8.1. PDAT.xxxx 

Data set PDAT.xxxx contains ply-level damage-dependent stresses, strains, and matrix crack internal state vari- 
ables. Data are centroidal values. The variable xxxx is the element name. The data for each element are stored in a 
record named DAT_PLY . ielt, where ielt is the element number. Each record contains n items, where 

n = nlayer * 9 

and nlayer is the number of layers in the model. 

The data are expressed with respect to ply coordinates and are stored in the following order: 

1 . Oi i is normal stress in the fiber direction. 

2. 022 > s normal stress transverse to the fibers. 

3. Oi2 is shear stress. 
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4. Eii is strain in the fiber direction. 

5. £22 is strain transverse to the fibers. 

6. E\2 is shearing strain. 

7. is internal state variable associated with fiber fracture. 

8. <*£22 * s i nterna l state variable associated with mode I opening of the matrix crack. 

9. <*£12 is internal state variable associated with mode II opening of the matrix crack. 

Repeated n layer times. 

B8.2. DGP.DATA.L1 

This data set is created by AUS/TABLE and contains the growth law parameters for the matrix cracking evolution- 
ary relationship. The following variables are used to specify table size: 

NI = number of material parameters, for this case 3 

NJ = number of material systems, for this case 1 

Type = numerical format, such as real or integer 

where NI and NJ are the number of columns and rows, respectively and Type specifies numerical format, real or integer. 
Each entry contains the following: 

1. Growth law parameter k. 

2. Growth law parameter h. 

d a L 

3. Parameter for determining -, dpara . 

dS 

These entries are repeated NJ times. 

B8J. ISV.xxxx 

This data set contains the matrix cracking internal state variables at each layer. The variable xxxx is the element 
name. The data are stored in a record named ALP AM. 1 . 

This record contains n items, where 

n = nlayer * nintgpt * nelt 

and nlayer is the number of layers in the model, nintgpt is the number of integration points for element, and nelt 
is the number of elements. 

The data are stored in the following order: 

1. is the internal state variable associated with fiber fracture. 

2. a L22 ls i nterna l state variable associated with mode I opening of the matrix crack. 

3. a^n is the internal state variable associated with mode II opening of the matrix crack. 

The data storage occurs for every layer, every integration point, and every element. 
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Appendix C 

Residual Strength Program 
Cl* General Description 

This appendix lists a sample program that was used to calculate the residual strength of a cross-ply laminate that was 
first fatigue loaded and then monotonically loaded to failure. The program is similar to that described in appendix B for 
Processor DGI. 

C2« Residual Strength Analysis Input 

The following list illustrates the input from a residual strength analysis. The problem being solved is the uniaxially 
tensile-loaded open-hole cross-ply laminated plate, which is shown in figure 6, and described in the main body of this 
report. The list contains the main program plus two procedure files. The first procedure file performs the calculations for 
each fatigue load cycle as described in appendix B. The second procedure file calculates the response during the mono- 
tonic loading to failure and is presented in this appendix. The finite element model was created using PATRAN. The file 
PT2T.CEHQUADFM.R1 .PRC was created with the PATRAN-to-testbed (PT2T) neutral file converter located in 
COMET. This file contains all the nodal locations, connectivity matrix, boundary conditions, and applied forces. 

#@$-o msg.out .Send output messages to file msg . out 

#@$-e msg. err .Send error messages to file msg. err 

#@$-q verylong@blackb . Batchfile queue 

#@$-me .Send mail when run is complete 

# 

cp $CSM_PRC/proclib . gal proclib.gal 

chmod u+w proclib.gal 

testbed > notchn.o << \endinput 

*set echo=of f 

*set plib=28 

*open 28 proclib.gal /old 
*open/new 1 cehquadatm. 101 

rectangular panel with circular cutout 
quarter panel mesh 
552 elements 
615 nodes 

EX4 7 4 node quad elements 

residual strength after fatigue and monotonic loading 
using monotonic growth law 

*add pf fb. clp 
*add pffdm.clp 
*def/a es_name = EX47 
*def/a es_proc = ESI 

*call ES ( function = 'DEFINE ELEMENTS'; es_proc = <es_proc> 
es_name=<es_name> ) 

[ xqt TAB 

START 615 

* ADD PT2T . CEHQUADFM . R1 . PRC .Runstream data from PATRAN modelling 
[xqt TAB 

*call PT2T_START .615 nodes 

j loc 

*call PT2T_JLOC .Obtain joint locations from PT2T.*.*.PRC 
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CONSTRAINT DEFINITION 1 
*cal 1 PT2T_BC 


. Constraints : 

.Fixed end and suppressed drilling dof from 
. PT2T . * . * . PRC 


[xqt AUS .Create input datasets 

SYSVEC : appl fore .Applied Forces 

*call PT2T_AF .Obtain applied forces from PT2T.*.*.PRC 

TABLE (NI=16,NJ=1) : OMB DATA 1 1 . Ply-level material property 
IM7/5260 
1=1,2, 3,4, 5 

J=l: 22.162E+6 0.333 1.262E+6 0.754E+6 0.754E+6 
1 = 6 , 7 , 8, 9 

J=l: 0.754E+6 1.0E-4 1.0E-4 0.01 
1=10,11,12,13,14,15,16 

J=l: 0.0 0.0 0.0 0.015 0.008 0.000 0.0 

TABLE (NI=3 , NJ=3 , i type=0 ) : LAM OMB 1 1 .Section properties 
J=1 : 1 0.006 0.0 
J=2 : 1 0.036 90.0 
J=3 : 1 0.006 0.0 


TABLE (NI=3 ,NJ=1, ITYPE=0) : DGP DATA 1 1 
IM7/5260 

J=1 : 1.1695 5.5109 3.8686E-7 

[xqt LAU 

ONLINE =2 
[xqt ELD 

*call PT2T_CONN 


[xqt E 

stop 

*call ES (function= ’ INITIALIZE 1 ) 

*call ES (function^ ’FORM STIFFNESS /MATL ' ) 
[xqt RSEQ 

reset maxcon=12 
[xqt TOPO 

reset maxsub=2 0 0000 
reset lram=100000 
reset lrkm=200000 
[xqt K 
[xqt INV 

*def/i ns__overwri t e=<t rue> 


.Damage evolution data 


.Create constitutive matrix 

.Define connectivity 
.Obtain connectivity from 
PT2T. * . * . PRC 

.Initialize initial datasets 

.Initialize element matrices 
.Form stiffness matrices 
. Resequence 

.Create maps 


.Assemble global stiffness matrix 
.Invert global stiffness matrix 


Call procedure to perform calculations at each cycle 

*call PFFB ( es_proc=<es_proc> ; es_name=<es_name> ; -- 

N_fcycl=l ; N__lcycl = 100 000 ; N_cylinc=20 .Fatigue up to 100,000 cycles 

by 20 cycle 
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NSUB=1 ; NSTRT=1 ; NS_lcycl = 50 ; -- .increments; ramp up in 

50 subincrements; 

NPRT=1000 ) .print datasets every 

1000^ increment 

Call procedure to perform monotonic loading 

*call PFFDM ( es_proc=<es_proc> ; es__name=<es_name> ; 

N_f eye 1-1 ; N_lcycl=2700 ; N_cylinc=l Increase load in 2700 

load steps by 1 step 

NSUB=0 ; NSTRT=0 ; NS_lcycl=0 ; -- .increment; no subincrements 

NPRT=100 ) .print datasets every 

lOO^ 1 increment 


*pack 1 
[xqt exit 
\endinput 

C2.1. Procedure to perform loop through calculations for each fatigue load cycle (file name pffb.clp) 

♦procedure PFFB { es_proc ; es^name ; 

N_f cycl ; N_lcycl ; N_cylinc ; 

NSUB ; NSTRT ; NS_lcycl ; NPRT ) 

Original version with subincrements 
Single major loop 

N_fcycl: first fatigue cycle 
N_lcycl: last fatigue cycle 
N__cylinc: cycle increment 
NSUB: subincrement flag 
NSTRT: cycle to start subincrements 
NS_lcycl : number of subincrements 
NPRT: output storage cycle increment 

begin loop here 

♦set echo=on,ma 
*set echo=off 
*def icount = 0 

*DO : CYCLOOP $NCYL = <[N_fcycl]>, <[N_lcycl]>, <[N_cylinc]> 

*def icount = ( <icount> + 1 ) 

*if < <icount> /eg <[NPRT]> > /then 
*def iprint = 1 
*def icount = 0 
*else 

*def iprint = 0 
*endi f 

*def $SNCYL = 1 

* IF < < < [NSUB] > /EQ 1> /AND < <$NCYL> /EQ < [NSTRT] > > > /THEN 
*def iscount = 0 
*DO $SNCYL = 1, <[NS_lcycl]> 

*def iscount = ( <iscount> + 1 ) 

*if < <iscount> /eq <[NPRT]> > /then 



*def isprint = 1 
*def iscount = 0 
‘else 

*def isprint = 0 
*endi f 

*def delinc = { 1.0 / <[NS_lcycl]> ) 

[xqt DRF 

select /PRINT = 0 
select /DSTATUS = 1 
select /XFACTOR = 0.0 
stop 

[xqt SSOL 

*if < <IPRINT> /eq 1 > /then 
[xqt VPRT 

format = 4 
print STAT DISP 
stop 

*endi f 

*call STRESS (directional; -- 

location= INTEG_PTS; print=<f alse> ) 

[xqt DGI 

select /PRINT = 0 

select /INC_SIZE = <delinc> 

select /N_CYCLE = <$NCYL> 

select /NINCR = <$SNCYL> 

select /DSTATUS = 1 

stop 

* i f < <ISPRINT> /eq 1 > /then 

*COpy 1, PLYDT . <ES_NAME> . <$NCYL> . <$SNCYL> 
1, PDAT . <ES_NAME> . * 

*COpy 1, DISP . <ES_NAME> . <$NCYL> . <$SNCYL> = 
1, STAT. DISP.* 

*copy 1, TISV . <ES_NAME> . <$NCYL> . <$SNCYL> 
1, ISV. <ES_NAME> . * 

♦copy 1, TSTRS . <ES_NAME> . <$NCYL> . <$SNCYL> 
1, STRS . <ES_NAME> . * 

*endi f 
*ENDDO 


♦JUMP : CYCLOOP 


♦ENDIF 




♦def 

delinc 

= < [N_cylinc ] > 

[xqt 

DRF 





select 

/PRINT = 

0 



select 

/DSTATUS 

= 

1 


select 

/XFACTOR 


0.0 


stop 




[xqt 

SSOL 




♦if < 

: <IPRINT> /eq 1 

> 

/then 



[xqt VPRT 

format = 4 
print STAT DISP 
stop 

*endi f 

*call STRESS (directional; location= INTEG_PTS; print =<false> ) 
[xqt DGI 

select /PRINT = 0 

select /INC_SIZE = <delinc> 

select /N_CYCLE = <$NCYL> 

select /NINCR = <$SNCYL> 

select /DSTATUS = 1 

stop 

*if < <IPRINT> /eq 1 > /then 

*Copy 1 # PLYDT . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, PDAT . <ES_NAME> . * 
*copy 1, DISP . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, STAT. DISP.* 

*copy 1, TISV . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, ISV . <ES_NAME> . * 
*copy 1, TSTRS . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, STRS . <ES_NAME> . * 

*endi f 

: CYCLOOP 

*set echo^off 
*end 

C2.2. Procedure to perforin monotonic loading calculations (file name pffdm.clp) 

^procedure PFFDM { es_proc ; es_name ; -- 
N_fcycl ; N_lcycl ; N_cylinc ; 

NSUB ; NSTRT ; NS_lcycl ? NPRT ) 

File to control monotonic loading to failure 
Original version with subincrements 
Single major loop 

N_fcycl: first load step 
. N_lcycl: last load step 

N_cylinc: load step increment 

NSUB: subincrement flag (=0, to bypass) 

NSTRT: step to start subincrement s (=0 , to bypass) 
NS_lcycl: number of subincrement s ( =1 , to bypass) 

NPRT: output storage step increment 

begin loop here 

*set echo=on,ma 
*set echo=off 
*def icount = 0 

*DO : CYCLOOP $NCYL = <[N_fcycl]>, <[N_lcycl]>, <[N_cylinc]> 

*def icount = ( <icount> + 1 ) 

*if < <icount> /eq < [NPRT] > > /then 
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*def iprint = 1 
*def icount = 0 
*else 

*def iprint = 0 
*endi f 

*def $SNCYL = 1 


*IF < < < [NSUB] > /EQ 1> /AND < <$NCYL> /EQ < [NSTRT] > > > /THEN 
*def iscount = 0 
*DO $SNCYL = 1 # <[NS_lcycl]> 

*def iscount = ( <iscount> + 1 ) 

*if < <iscount> /eq <[NPRT]> > /then 
*def isprint = 1 
*def iscount = 0 
*else 

*def isprint = 0 
*endif 

*def delinc = ( 1.0 / <[NS_lcycl]> ) 

[xqt DRF 

select /PRINT = 0 
stop 

[xqt SSOL 

*if < <IPRINT> /eq 1 > /then 
[xqt VPRT 

format = 4 
print STAT DISP 
stop 

*endif 

*call STRESS ( direct ion=l ; -- 

location^ INTEG_PTS; print=<f alse> ) 

[xqt DGI 

select /PRINT = 0 
select /INC_SIZE = <delinc> 
select /N_CYCLE = <$NCYL> 
select /NINCR = <$SNCYL> 
stop 

*if < <ISPRINT> /eq 1 > /then 

*copy 1, PLYDTM . <ES_NAME> . <$NCYL> . <$SNCYL> = -- 
1, PDAT . <ES_NAME> . * 

*copy 1, DISPM . <ES_NAME> . < $NCYL> . < $SNCYL> = -- 
1, STAT. DISP.* 

*copy 1, TSTRS . <ES_NAME> . <$NCYL> . <$SNCYL> = -- 
1, STRS . <ES_NAME> . * 


*endif 

*ENDDO 

* JUMP : CYCLOOP 
*ENDIF 



*def delinc = <[N_cylinc]> 

[ xqt DRF 

select /PRINT = 0 
select /DSTATUS = 22222 
select /XF ACTOR = 0.00079 
stop 

[xqt SSOL 

*if < <IPRINT> /eq 1 > /then 
[xqt VPRT 

format = 4 
print STAT DISP 
stop 

*endi f 

♦call STRESS (direction=l ; location= INTEG_PTS ; print=<false> ) 

[xqt DGI 

select /PRINT = 0 
select /INC_SIZE = <delinc> 
select /N_CYCLE = <$NCYL> 
select /NINCR = <$SNCYL> 
select /DSTATUS = 22222 
stop 

*if < <IPRINT> /eq 1 > /then 

♦copy 1, PLYDTM . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, PDAT . <ES_NAME> . * 

♦copy 1, DISPM . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, STAT. DISP.* 

♦copy 1, TISV. <ES_NAME> . <$NCYL> . <$SNCYL> = 1, ISV . <ES_NAME> . * 

♦copy 1, TSTRS . <ES_NAME> . <$NCYL> . <$SNCYL> = 1, STRS . <ES_NAME> . * 

♦endif 

: CYCLOOP 

♦set echo=off 
♦end 

C3. Residual Strength Analysis Output 

The following lists illustrate the standard output from a residual strength analysis. The print flag is set equal to 0 so 
that the only information stored in the output file is the cycle number, failed ply number, the current eji, and the current 
Oj i for the failed elements. The stress, strain, and displacement fields are still stored in the library data sets as are the 
internal state variables. How often such data are stored in data sets is up to the user and is controlled by the NPRT vari- 
able in the runstream and the *copy 1 command in the procedures pffb.clp and pf fdm. clp. 

The first list is at the end of the fatigue loading, cycle number 99981 . 


** BEGIN DGI ** DATA SPACE= 2000000 WORDS 


CYCLE 
+ + PLY 

NUM. 
1 OF 

ELEMENT 

99981 

8 

0 . 1832E-01 

0 . 2743E+06 

99981++ 

+ + PLY 

3 

OF 

ELEMENT 

8 

0 . 1832E-01 

0 . 2743E+0 6 

9 9 9 81++ 

+ + PLY 

1 

OF 

ELEMENT 

9 

0 .2511E-01 

0 . 3 023E+0 6 

9 99 81++ 

+ + PLY 

3 

OF 

ELEMENT 

9 

0 .2511E-01 

0 . 3 023E+0 6 

99981++ 
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++PLY 1 OF ELEMENT 
++PLY 3 OF ELEMENT 
++PLY 1 OF ELEMENT 
++PLY 3 OF ELEMENT 
++PLY 1 OF ELEMENT 
++PLY 3 OF ELEMENT 
++PLY 1 OF ELEMENT 
++PLY 3 OF ELEMENT 
EXIT DGI CPUTIME= 

CONVEX COMET VER . 1 . ! 


10 0.4472E-01 

10 0 . 4472E-01 

18 0 . 1554E-01 

18 0 . 1 554E-01 

19 0 . 1817E-01 

19 0 . 1817E-01 

20 0 . 1948E-01 

20 0 . 1948E-01 

1.74 

.4 - DEC. 1994 


0 . 3324E+ 06 

99981++ 

0 . 3324E+06 

99981++ 

0 . 3242E+06 

99981++ 

0 . 3242E+06 

99981++ 

0 . 3324E+06 

99981++ 

0 . 3324E+06 

99981++ 

0 . 3324E+06 

99981++ 

0 . 3324E+06 

99981++ 


(blackb) 


07:19:95 18:58:26 


The second list is at load step 715 in the monotonic loading procedure. Since xf actor = 0.00079 in pffdm.clp, and the 
applied load is 1572 lb/in, the load step 715 corresponds to applied load + x factor * applied load *load step = 2460 
lb/in. 


** BEGIN DGI ** DATA SPACE= 2000000 WORDS 


CYCLE 

NUM. 

= 

715 




+ + PLY 

1 

OF 

ELEMENT 

7 

0 . 1542E-01 

0 . 2988E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

7 

0 . 1542E-01 

0 . 29 88E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

8 

0 .1965E-01 

0 . 3040E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

8 

0 . 1965E-01 

0 . 3040E+06 

7 1 5 + + 

+ + PLY 

1 

OF 

ELEMENT 

9 

0.253 0E-01 

0 . 3067E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

9 

0 .2530E-01 

0 . 3 067E+ 06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

10 

0 .2093E+00 

0 . 3382E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

10 

0 .2093E + 00 

0 . 3382E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

18 

0.1674E-01 

0.3209E+06 

715++ 

+ + PLY 

3 

OF 

ELEMENT 

18 

0 . 1674E-01 

0.3209E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

19 

0.3489E-01 

0 . 32 90E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

19 

0 .3489E-01 

0 . 32 9 0E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

20 

0 . 13 1 6E+00 

0 . 33 84E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

20 

0.1316E+00 

0 . 33 84E+06 

7 15 + + 

+ + PLY 

1 

OF 

ELEMENT 

28 

0.172 8E- 0 1 

0.3279E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

28 

0.1728E-01 

0 . 3279E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

29 

0 .3099E-01 

0.3297E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

29 

0 . 3 099E-01 

0 . 32 97E+ 06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

30 

0.8494E-01 

0 . 33 89E+06 

715++ 

+ + PLY 

3 

OF 

ELEMENT 

30 

0 . 8494E-01 

0.3389E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

38 

0 . 1595E-01 

0.3291E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

38 

0 . 1595E-01 

0 . 32 91E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

39 

0 .2312E-01 

0 . 32 97E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

39 

0.2312E-01 

0 . 32 97E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

40 

0.5433E-01 

0 . 33 66E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

40 

0 . 5433E-01 

0.33 66E+06 

7 15 + + 

+ + PLY 

1 

OF 

ELEMENT 

48 

0 . 1555E- 01 

0 . 33 17E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

48 

0.1555E-01 

0.3317E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

49 

0 .2130E-01 

0 . 32 92E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

49 

0 .2130E-01 

0.3292E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

50 

0 .2729E-01 

0.3337E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

50 

0.2729E-01 

0 . 3337E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

59 

0 . 1883E-01 

0.3333E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

59 

0 . 1883E-01 

0.3333E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

60 

0.1919E-01 

0 . 332 9E+06 

715 + + 

+ + PLY 

3 

OF 

ELEMENT 

60 

0 .1919E-01 

0 . 3329E+06 

715 + + 

+ + PLY 

1 

OF 

ELEMENT 

70 

0 . 169 9E-01 

0 . 33 12E+06 

715 + + 
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++PLY 3 OF ELEMENT 70 0 . 1699E-01 0.3312E+06 715++ 

EXIT DGI CPUTIME= 1.68 

CONVEX COMET VER. 1.5.4 - DEC. 1994 (blackb) 07:19:95 22:23:34 
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Table 1. Material Properties of Unidirectional Ply of IM7/5260 


£ u ,Msi 22.16 

£ 2 2 . Msi 1.26 

Gj 2 , Msi 0.75 

Vj2 0.333 

/pl y , in 0.006 

Silent °' 015 

S 22 cnt 0008 

Growth law parameters: 

k 1.1695 

n 5.5109 

dpara 3.8686x 10 


Table 2. Maximum Fatigue Loads Employed in Sample Calculations 


Layup 

Specimen geometry 

Maximum fatigue 
load (/? = 0.1), Ib/in. 

[0/±45/90] s 

Unnotched 

3300 

Open hole 

2000 

[0/90 3 ] s 

Unnotched 

2480 

Open hole 

1572 





Figure 1. Progressive failure analysis scheme. (From ref. 13.) 
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Figure 5. Global displacements from load redistribution. 


Nx N x 



Loading direction 



y 


Figure 7. Finite element model for a laminate with a central circular hole. 
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1.02 


Normalized 
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Figure 8. Stiffness loss of IM7/5260 laminates with central circular hole. 
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6000 
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Figure 9. Predictions of residual strength. 
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til ‘ i 1 i i i i lit ii lii it lii ii lii ii I i i n,l i i ii I 

0 .05 .10 .15 .20 .25 .30 .35 .40 

Distance from notch, in. 

(a) No damage. 



Distance from notch, in. 

(b) Ply discount, y = 0; Ultimate tensile load = 2379 lb. 



Distance from notch, in. 


(c) Elastic perfectly plastic, y = 1 ; Ultimate tensile load = 3435 lb. 
Figure 1 0. Fiber failure criteria. 
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(b) Mesh 2. 1071 nodes; 992 elements. 



(c) Mesh 3. 2225 nodes; 2112 elements. 
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(d) Mesh 4. 2813 nodes; 2688 elements. 

Figure 1 1 . Finite element meshes used in convergence study. 


40 


Ultimate 
ten: 
loac 


387 nodes 1071 nodes 2225 nodes 2813 nodes 

336 elements 992 elements 2112 elements 2688 elements 

Figure 12. Mesh refinement study for residual strength predictions of [0/±45/90] s laminate open-hold geometry 
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